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We generalize Greenberger-Horne-Zeilinger (GHZ) theorem to an arbitrary number of D- 
dimensional systems. Contrary to conventional approaches, we use concurrent observables, which 
are incompatible but still have a common eigenstate. We begin with the theorem for 4 systems of 
a dimension divisible by 3, and discuss its extension to TV systems of an arbitrary dimension. The 
GHZ theorem can be proved as long as TV is not divisible by all non-unit divisors of D, smaller than 
TV. 



I. INTRODUCTION 

Bell's theorem says there exists quantum correlations 
between separated observers which cannot be simulated 
by local realistic theories Using a statistical in- 

equality, Bell showed that any local realistic theories are 
incompatible with quantitative predictions of quantum 
mechanics for a singlet state of two spin- 1/2 particles [1[. 
Even though Bell's theorem was studied mostly in terms 
of statistical inequalities, different version of it without 
inequalities, was also shown for a multiqubit systems by 
Greenberger, Home, and Zeilinger (GHZ) Q. They de- 
rived an all-versus-nothing contradiction based on per- 
fect correlations of the so-called GHZ states. This leads 
to a direct refutation of EPR ideas on relation between 
locality and elements of reality with quantum mechanics. 

Extending Bell's theorem is important not only for 
a deeper understanding of foundations of quantum me- 
chanics. It allows developing new applications in quan- 
tum information processing, such as quantum cryptog- 
raphy, secret sharing, quantum teleportation, reduction 
of communication complexity, quantum key distribution 
and random numbers generation [5|-l9|. A Bell- type test 
is also regarded as a fundamental method to verify en- 
tanglement among subsystems. Therefore, search for ex- 
tensions of Bell's theorem to more complex systems such 
as multipartite and/or high-dimensional systems is an 
important task in quantum information science [3|- 

An all-versus-nothing test, that we call GHZ theorem, 
has been generalized to higher dimensional systems. For 
the sake of convenience, we shall use the tuple (TV, M, D) 
to denote TV parties, M measurements for each party and 
D distinct outcomes for each measurement. In Ref. fioj ] . 
GHZ theorem was derived for a (D + 1, 2, D) problem. A 
probabilistic but conclusive GHZ-like test was shown for 
(D,2,D) in Ref. [11]. The (TV, M, D) problem for odd 
TV > D, M = 2, and even D was studied using operator 
relations by Cerf et al. [l2|. Recently, Lee et al. em- 
ployed incompatible composite observables for arbitrary 
(odd TV, 2, even D) problems [13|. The conventional ap- 
proaches are based on compatible observables. 

We generalize the GHZ theorem to 3 or more D- 



dimensional systems. To this end, we employ concurrent 
composite observables. Such observables are mutually 
incompatible but still have a common eigenstate, here 
a generalized GHZ state. They can be realized by using 
multiport beam splitters and phase shifters, as it is shown 
in Refs. [13, H, G3 • We illustrate our principal idea with 
four 3d-dimensional systems, and provide a systematic 
method to extend it to three or more D-dimensional sys- 
tems. As non-unit divisors of D play a crucial role, one 
can show a GHZ-type contradiction as long as TV is not 
divided by all non-unit divisors of D, which are smaller 
than TV. We also show that our approach can reproduce 
the previous results 0, OH EH Gil an d it can lead to a 
general GHZ theorem for TV quDits. 



II. CONCURRENT OBSERVABLES AND 
INVARIANCE CONDITION 

Some sets of observables have a common eigenstate. 
If a system is prepared in the eigenstate, the measure- 
ment results for such observables are concurrently ap- 
pearing with certainty. Such observables are called "con- 
current" [13]. For a quantum system of dimension D, 
greater than two, consider two Hermitian operators A 
and B, such that A = a\ijj){^\ +A' and B = b\^){^\ +B'. 
The state is then a common eigenstate of both observ- 
ables as A{B)\ijj) = a(b)\^), even if [A,B] = [A', B'\ ^ 0. 
Note that compatible observables are clearly concurrent. 

To construct concurrent observables, in a Hilbert space 
ft, one can use the method of Ref. [13|- Consider 
a unitary operator £7, which is of the form of U = 
e^l^^l + C/ 1 - with £7^) =0. Here U 1 - is a unitary op- 
erator on a space 7~L ± which is denned by the requirement 
n = n^ ©ft- 1 , where is the one-dimensional space 
containing Every such unitary operator leaves the 
state unchanged, up to a global phase. If the state 
|^) satisfies A\ijj) = A|^), then all transformed operators 
B\j = UAU^ are concurrent with A. 

We shall use the following local observables. Con- 
sider a measurement in .D-dimensional Hilbert space 
represented by A = ^2^0 ujn \ n ) AA ( n \^ wnere 00 = 



exp(27rz / D) . The operator A is unitary, however one can 
uniquely relate it with a Hermitian observable H, by re- 
quiring A = exp(iH). Therefore the complex eigenvalues 
of A, which could be called a unitary observable, can be 
associated with the measurement results, denoted by real 
eigenvalues of H [HI, [l3|, IH| . Such a unitary representa- 
tion leads to a simplification of the mathematics without 
changing any physical results. 

Consider N quDits in a generalized GHZ state, 



J 2 (n) 



D-l 



® |n), 



(1) 



n=0 



TV 



where \n) denotes a basis state for a quDit. Let us con- 
sider a unitary operation in the form ofU = (S)k=i -^k(fk) 
with phase shifter P k = Y^n=o \ n )( n \- ^ "phases" 

fk(n) satisfy a condition 



TV 



mod D, 



(2) 



k=l 



for each n, then the unitary operator U leaves the GHZ 
state |^) invariant. 

Consider now a composite observable X® N = X ® 
X <g> • — <g> X, where the local observable X is defined 
by applying quantum Fourier transformation F on a ref- 
erence unitary observable Z = Yln=o ^ n |^)(^|? that is 
X = FZF^ . An eigenvector of X associated with the 
eigenvalue uo n is given by 



n 



F\n) 
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(3) 



With the standard basis set {|n)}, the observable X can 

put as by X = Yln=o \ n )( n + Mi where \n) = \n mod D). 
The generalized GHZ state is an eigenstate of the 
composite observable X® N with the unity eigenvalue: 
X® N \ip) = X <g> X <8> • • • <8> i#) = 1^). 

We can construct a set of concurrent composite observ- 
ables by applying the unitary operation U = (S)fcLi Pk 
to X®^. Consider the phases /^(n) = a^n with ratio- 
nal numbers a k . If the phases fk(n) satisfy the invari- 
ance condition (|2]), the transformed composite observable 
(jX®N(j] = x(ai) 0X(a 2 ) • • >®X(a N ) is concurrent 
with i.e., UX® N U*\i)) = 

For each eigenvalue o; n , the eigenvector of the local 
observable X(a) is given by applying the phase shifter P 
on \n) x as 



\n) a P\n) x 



1 



D-l 



(n+o;)m 



m). 



(4) 



The observable X(a) can be written in the standard basis 
set {\n}} as 



m-2 



X{a) = u,-* [J2 Mi* + M+ " aD \D " 1>(0| • (5) 
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FIG. 1. The invariant condition for 4-partite system. Each 
party applies a local phase shifter Pk with phase fk (n) to the 
generalized GHZ state \ip). If the sum of the local phases 
is congruent to zero modulo D, i.e., fi(n) + /2(n) + 
h(n) + f±(n) = mod then the state transformed by 
the phase shifters is unchanged. 



Note that if a is an integer, the measurement basis set 
of X(a) will be the same as {|n) x } of X except the 
ordering, i.e. \n) a = \n + a) x , and then X(a) = uo~ a X. 
That is, the observable X(a) is equivalent to X, up to 
a phase factor uj~ a . Let \n) a be the eigenstate of X(a) 
associated with eigenvalue cj u , and \m)p the eigenstate 
of X(f3) associated with uj m . If and only if a differs 
from f3 by an integer, then two measurement bases sat- 
isfy \ a (n\m)f3\ 2 = where 7 = m — n + /3 — a. Here 
^d(t) = 1 if 7 is congruent to zero modulo D and other- 
wise Sd(j) = 0. That is, unless f3 — a is an integer, two 
local observables X(a) and X(f3) are inequivalent. 



III. GENERALIZED GHZ THEOREM 



A. Four-partite systems 



We illustrate our idea by considering first a four qu£ut 
system. Already this case goes significantly beyond the 
previous studies [To|, EH Elf. 

Take a 4 quDit GHZ state (D is an integral multiple 
of 3) 



1 D_1 

W) -7= Yl \ n ^ n ^ n ^ n ) 

V D n =0 



(6) 
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The quDits are distributed to four sufficiently separated 
parties. Each party performs one of two non-degenerate 
local measurements on his or her quDit, each of which 
produces distinguishable D outcomes. The eigenvalues 
of the observables are of the form uj k . 

Consider such a measurement, Mi for the i-th 
party. One can denote the probability that the par- 
ties obtain results a; mi , where rrii is an integer, by 
V(mi, ?7i2, 7713, 7714), and define a correlation function 
E QM (M 1 ,M 2 ,M 3 ,M 4 ) = {^\M X ® M 2 ® M 3 ® M 4 |^>, 
the quantum average of products of the measurement re- 
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suits: 



D-l D-l 

"' J2 ^ mi+m2+m3+m4 ^(^i,-^4). (7) 

mi=0 772,4=0 



If local realism is also to reproduce quantum perfect 
correlations for the case of vq = X® 4 , one must have 



Elr(xi,X2,x 3 ,X4 : ) = 1, 



If all measurements are X, that is each M ? = X, since or f° r eacn run 



X®X®X®X\tl>) = |^>, 



(8) 



one has £?qm(X, X, X, X) = 1. We have a perfect GHZ 
correlation. If arbitrary three parties know their own 
outcomes u Xi of measurements X, then they can predict 
with certainty the remaining party's outcome: We will 
denote such a perfect correlation by 



Cqm(%1 + X 2 + Xs + X4 = 0). 



(9) 



Of course, the sum is taken modulo D; such a convention 
is used below in all formulas. 

We can now construct concurrent composite observ- 
ables from the observable % = X® 4 , by applying a uni- 
tary operator Ui = A ® P 2 ® 3 , w ^ n the phase shifters A 
of phases fi(n) = (D — l)n and /2(w) = n/3. One of the 
new observables is «i = E/iiWi = X(£> - 1) (g)X(l/3) 03 . 
The phases //c(n) satisfy the invariance condition (|2]), 
fi(n) +3/2(n) = mod D. Thus, the observable #i has 
1*0) as its eigenstate with eigenvalue 1. By Eq. ([5]) one 
has X(D - 1) = cjX(O), and Y = X(l/3) is given by 



m-2 



Y 



-1/3 



^|n)(n + l|+^/ 3 ^-l)(0| . (10) 
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The observable Y = X(l/3) is not equivalent to X = 
X(0), see the explanation below Eq. ((5]). 

The other three concurrent observables are obtained by 



UiVqUJ , Z G {2,3,4}, where the unitary operators 
Ui are composed of the cyclic permutations of the local 
phase shifters: U 2 = A® A® A® A, Us = A® A® A® 
A ? an d C/4 = A ® A ® A ® A- The perfect correlations 
for ^ (i = 1, . . . , 4) are respectively given by 

C QM (xi + y 2 + 2/3 + 1/4 = -1), 
Cqm(2/i + ^2 + 2/3 + 1/4 = -1), 
Cqm(vi + 2/2 + ^3 + 2/4 = -1), 
Cq M (2/i + V2 + 2/3 + ^4 = -1), (11) 

where uo Vi is an outcome of measurement Y for i-th party. 

Local realistic theories assume that the outcomes of the 
measurements are predetermined, before the actual mea- 
surements. This means that the values of local realistic 
predictions for the correlations ([TT]h for each experimen- 
tal run, must satisfy: 



(J yi UJ y2 UJ y3 UJ X4 = 



UJ 



UJ 



UJ Xl UJ X2 UJ X3 UJ X4 



1, 



(13) 



(14) 



(12) 



If one multiplies side-by-side all equations (fT2j) . one sees 
that this is possible only provided oj~ 3 ^^ yi ~ 4z = 1. As 
uo = exp(27ri/D), if D = 3d where d is an integer, an ele- 
mentary algebra shows that there is no integer solution of 
y of the equation 3y + 4 = mod D. Thus local realistic 
correlation (p3J) is impossible (a GHZ-like contradiction). 

It is worth noting that the approach with concurrent 
observables relaxes the restrictions of early studies requir- 
ing compatible observables. This enables one to generalize 
GHZ contradictions, beyond the case N > D studied in 
Refs. HQ. 

Note that in order to prove the 4-partite GHZ contra- 
diction we have chosen the local dimension D and the 
number of the observables Y's in the considered cor- 
relation functions, which will be denoted as X2, such 
that the greatest common divisor (gcd) of D and 7V 2 , 
gcd(Z},iV2) = 3, does not divide the number of parties 
N = 4, or equivalent ly the number N± of the observables 
X (here equal to 1). This mathematical property plays a 
central role in the generalization of the GHZ contradic- 
tion to an arbitrary number of parties. 

The unitary observables X and Y can be realized by 
using optical devices such as multiport beam splitters 
and phase shifters (see, [l6[). 



B. Extension to multipartite systems 

We extend our approach a general case of N quDits, 
which N > 3, which is non-divisible by any non-unit 
divisor of D, smaller than N. 

To this end, we use a set of (N + 1) concurrent ob- 
servables given by £0 = X® N and N observables of 
the following forms: v x = uuX® Nl <g> Y® N2 and v k = 
y®fc-i ^x^n, ^ y-0A/- 2 -/c+i for k = 2 ,...,X 2 + 1 
and finally v k = X^-n 2 -i f®N 2 ^x^N-k+i for 
/e = X 2 + 2,...,X 

The composite observable V\ is obtained by a unitary 
transformation, U ± = A ® 1 0A/W <g> P£> N2 , o f the ob- 
servable vq, i.e., v\ = U\VoU\, with the phase shifters A 
and A of phases fi(n) = (D — l)n and f2(n) = n/7V2, 
respectively. The local observables X = X(0) and 
Y = X(l/N2) are given by Eq. ([5]). Likewise, we obtain 
the other concurrent observables vi (2 < I < N) by cyclic 
permutations in the unitary operator ZTi, as it was done 
in Sec. IIIII The phases satisfy the invariance condition in 
Eq. (|2J) as /i(n) + X 2 / 2 (n) = mod D. Thus, the TV- 
partite generalized GHZ state = -i= Yln=o \ n )® N 1S 



a common eigenstate of all the (7V + 1) concurrent observ- 
ables vi (1 = 0,..., N), with the same eigenvalue 1. This 
leads to the following values of correlation functions (for 
later convenience we now give to the parties the indices 
i = 1, ...TV, which will later on allow us to get a more con- 
cise notation in formulas). For all local observables equal 
to X, that is global Vq, one has Eqm (X, X, . . . , X) = 1. 
Thus we have a perfect correlation which can be denoted, 
in the way introduced earlier, as Cqm(X^i x i =0)- For 
Vk, where k = 1,2,. ..TV, one has perfect correlations of 
the following forms: for k = 1 

iVi N 

Cqm(X^+ ^ = -!)> ( 15 ) 

2=1 i=N 1 + l 

and for k = 2,...,7V 2 + 1 

k-1 iVi+fe-l N 

cqm£^+ Yl Xi + ^ = -1 )' ( l6 ) 

2=1 i=/c i=N!+k 

and finally for fc = 7V 2 + 2, iV 

k-N 2 -l k-1 N 

2=1 i=k—N2 i=k 

Following the similar arguments as in the case of 4- 
partite GHZ contradiction, we obtain the following con- 
dition for the local realistic correlation function as for the 
composite observable % = X® N , to have value equal to 
the quantum prediction, that is 1. It reads (modulo D) 

N N 

2=1 2=1 

However, if 7V 2 is an integral multiple of g but N is not, 
i.e., N2 = rjg, and TV cannot be divided by g, then there 
are no solutions of y = yi to the equation iV 2 y+7V = 
mod D. The greatest common divisor of 7V 2 and D is 
an integral multiple of g, i.e., gcd(7V 2 , D)=kg for some 
positive integer k but kg cannot divide TV as N is not an 
integer multiple of g. Thus we have a contradiction with 
the quantum prediction. 

In summary, in order to show a GHZ contradiction for 
TV-partite and D-dimensional system, we choose that (a) 
D = dg, (b) iV 2 = rjg, where d and 77 are positive integers, 
and (c) N cannot be divided by g. 

The integer g, a non-unit divisor D, plays a crucial 
role. For example, let us consider four 6-dimensional 
systems. The non-unit divisors g, smaller than N = 4, 
are 2 and 3. If we choose g = 2, then we are unable to see 
any 4-partite GHZ contradiction as the greatest common 
divisor (gcd) of N 2 and D, gcd(7V 2 = 2,D = 6) = 2, 
also divides N = 4. If we choose g = 3, on the other 
hand, gcd(7V 2 = 3, D = 6) = 3 but g = 3 does not 
divide N = 4 and the 4-partite GHZ contradiction can 
be proved. This is in fact a specific example of GHZ 



contradiction for a (4,2,3d) problem. In other words, 
we are always able to prove the GHZ contradiction for 
N(> 3)-partite and D(> 2)-dimensional system as long 
as N cannot be divided by all non-unit divisors of D. 

C. Comparison with previous works 

For appropriate values of iV and D , ou r approach re- 
produces the previous works [To|, EH, 13 . The work m 
Ref. [Toj| showed a GHZ contradiction for (D + l) qu^Dits. 
In our scheme, we are able to choose the number Ni = 1 
of the observable X and 7V 2 = D of the observable Y, 
then our GHZ contradiction for (N,2,D) system is re- 
duced to a (D + 1,2,D) problem; N± = 1 is indivisible 
by D = gcd(iV 2 = D,D). Consider a non-unit divi- 
sor g = 2 of D and an arbitrary odd integer N±. Thus 
we have a (odd 7V, 2 , even D) problem, which was stud- 
ied in the Refs. [12|, [l3|. One can also easily check that 
if 7V 2 = D = 2 and N± = 1, then our contradiction 
is reduced the original one of Greenberger, Home, and 
Zeilinger Q. 



IV. GENUINELY MULTIPARTITE 
HIGH-DIMENSION GHZ THEOREM 

The GHZ theorem for 2-dimensional systems seems to 
be fully understood. However, for more complex sys- 
tems this is not so. Cerf et ai suggested a criterion for 
a genuinely iV-partite (D-dimensional) GHZ contradic- 
tion [12]. It must arise only for the given TV-partite (D- 
dimensional) system, but not for the n(< 7V)-partite sub- 
set, or for effectively lower dimensionality of the involved 
observables (example: the three qubit classic GHZ theo- 
rem can put as theorem for three qutrits, and specific en- 
tangled GHZ states involving only two dimensional sub- 
spaces for each qutrit) . 

The GHZ contradiction we show here is a genuinely 
A/"-partite one, as it is constructed using a set of compos- 
ite observables composed of cyclic permutations. Let us 
explain this with the 4-partite GHZ contradiction. We 
use the following five concurrent composite observables: 

X <g> X <g> X <g> X 
X®Y®Y®Y 
Y®X<£)Y®Y 
Y (8) y <g>X ®Y 

Y®Y®Y®X (19) 

In such a circumstance, if we eliminate one of the parties, 
we are not able to show a GHZ contradiction with the 
remaining observables. The 4-partite GHZ state is no 
longer their common eigenstate. Similar argument can 
be put forward in the case of our iV-partner theorem. 

The genuine ^-dimensionality is reflected by the fact 
that the operators are undecomposable to a direct sum 
of any subdimensional observables [13|- In other words 
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if two local observables X and Y can be simultaneously 
block-diagonalized by some similarity transformation S 
such that SX& = X x © • • • © X K and SYS* = Yi © 
• • • © Yx-, then there exist some eigenstates \n) a of X and 
\m)p of Y such that a (n\S^ S\m) p = and one can find a 
sub-dimensional GHZ contradiction. However, there are 
no such eigenstates because for every n and m 



\a(n\m)p\ : 



sin 2 (7r£) 
D 2 sin 2 [(tt/D)^} 



>0, 



(20) 



where £ = m — n + /3 — a. As, the local observables 
X = X(a) and Y = X(j3) are such that /3 — a is not an 
integer, £ is a non-integral rational number. Thus, our 
GHZ contradiction is genuinely D-dimensional. 

One might try to generalize our method to more than 
two measurement settings. However this cannot be a 
genuinely multi-setting case. Consider a (TV, M, D) = 
(5,3,2) problem. Our approach requires six concurrent 
observables v = X©X©X©X©X, v ± = ujX <g>Y <g>Y © 
Z © Z, and vi formed by cyclic permutations of observ- 
ables in vi(2 < I < 5). The local observables are given by 
the Eq. (J5J) with the f^n) = (D - l)n, f 2 (n) = n/8, and 
H( n ) = 3n/8, respectively. As the phases satisfy the in- 
variance condition in Eq. (j2j) as /i(n) + 2/2(71) + 2/3(71) = 
mod Z}, the 5-partite and 2-dimensional GHZ state 
\ip) = (|0, 0, 0, 0, 0) + 11, 1, 1, 1, l))/<\/2 is a common eigen- 
state of these concurrent observables, with the eigenvalue 
equal to 1. With all that, one can show a GHZ contra- 
diction. However, it is not a genuine 3-setting one. The 
contradiction also arises for observables, composed of two 
local measurements X and Y, uq = X©X©X©X©X, 
ui = uoX ©l 7 '©! 7 '©! 7 "©^, and ui formed with cyclic 
permutations in u\(2 < I < 5), of a similar kind as above. 
Thus problem whether there are genuinely multi-setting 
GHZ contradictions is still open. 



V. SUMMARY 

We gave a generalized GHZ contradiction for multi- 
partite and high- dimensional systems. For the purpose, 
we employed concurrent composite observables, which 
have a generalized GHZ state as a common eigenstate 
(even though these observables are incompatible). We 
obtained concurrent observables by using the unitary op- 
erations composed of the local phase shifters. Previous 
approaches had used compatible composite observables. 
Abandoning this allowed to reach a broader class of GHZ 
contradictions. 

As a result, we showed a 4-partite GHZ contradiction 
for 3d-dimensional systems, where d is an arbitrary in- 
teger. We extended the method to Af-partite situation 
with D-dimensional systems. It turned out that non-unit 
divisors of D plays a crucial role to show the GHZ con- 
tradiction. We were always able to show a contradiction 
as long as N cannot be divided by all non-unit divisors 
of smaller than N. We also showed that our formu- 
lation of a generalized GHZ contradiction is genuinely 
Af-partite and D-dimensional. 
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